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We discuss the electronic structure, lattice dynamics and electron-phonon interaction of newly 
discovered superconductor LaOo.5Fo.5BiS2 using density functional based calculations. A strong 
Fermi surface nesting at k=(7r,7r,0) suggests a proximity to charge density wave instability and 
leads to imaginary harmonic phonons at this k point associated with in-plane displacements of 
S atoms. Total energy analysis resolves only a shallow double-well potential well preventing the 
appearance of static long-range order. Both harmonic and anharmonic contributions to electron- 
phonon coupling are evaluated and give a total coupling constant A ~ 0.85 prompting this material 
to be a conventional superconductor contrary to structurally similar FeAs materials. 

PACS numbers: 74.20.Pq, 74.70.-b 



Superconductors with layered crystal structures such 
as cuprates[l[, ruthenatesjH , or MgB2@ have gener- 
ated enormous research interest. A recent discovery 
of iron pnictidesjj] have triggered another wave of ex- 
tensive studies 5], and while the mediator of pairing 
in these systems remains officially unidentified, a large 
amount of evidence points to magnetic spin fluctuations 
induced by antiferromagnetic spin-density-wave (SDW) 
instability due to Fermi-surface nesting |6j at wave vec- 
tor k=(7r, 7r, 0) similar to the cuprates. Usually changing 
the blocking layer can tune the superconducting T c , thus 
searching for new layered superconductors is of both fun- 
damental and technological importance. 

Very recently, a new layered superconductor Bi404S3 
has been foundR, and soon after, two other systems, 
LaOi_ a; F :r BiS2 [Hand NdOBiS 2 [g] have been discov- 
ered. Here, the basic structural unit is the BiS2 
layer which is similar to the Cu-0 planes in Cu-based 
superconductors [l| and the Fe-A (A=P, As, Se, Te) 
planes in iron pnictidesj^]. A chance to explore supercon- 
ductivity and increase T c in these new compounds has 
already resulted in a lot of work that appeared shortly 
after the discovery P, @, 9 16]. Hall effect measurements 
reveal multiband features and suggest the superconduct- 
ing pairing occurs in one-dimensional chains [1 0t | . It was 
proposed that these compounds are type II supercon- 
ductors and good candidates for thermoelectric materi- 
als 
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Electrical resistivity measurements under pres- 
sure reveal that Bi404S3 and La(0,F)BiS2 have differ- 
ent T c versus pressure behavior, and the Fermi surface 
of La(0,F)BiS2 may be located in the vicinity of some 
instabilities 12] . A two p bands electronic model has been 
proposed based on band structure calculation fl3jj . and a 








FIG. 1: (a) Structure of La(Oo.sFo.5)BiS 2 ; (b) unstable 
phonon mode corresponding to wave vector k=(7r,7r,0). 



good Fermi-surface nesting with wave vector k=(7r, 7r, 0) 
has been found (l3|. The importance of the nesting has 
been emphasized experimentally [l2j and it was suggested 
that electronic correlations may play a role in supercon- 
ductivity of these systems [3] ■ 

Here we report our theoretical studies of the elec- 
tronic structure and lattice dynamic properties for 
LaO .5F .5BiS 2 , the compound that posses the highest 
T c ~10 K Q among known BiS based materials^, 13, % 
and whose structure is similar to superconducting iron 
arsenides LaFeOi_ x F K As 4]. Our first-principles calcu- 
lations are based on density functional theory (DFT) and 
show that the bands around the Fermi level are not sen- 
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TABLE I: Calculated lattice parameters and Wyckoff posi- 
tions of LaOBiS2 and La(Oo.5Fo.s)BiS2. Experimental results 
[l| are also listed for comparison. 
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FIG. 2: Band structure of (a) La0 .5Fo.5BiS 2 ; (b) LaOBiS 2 



sitive to F substitution, thus making rigid band approx- 
imation adequate for LaOi_ a; F a ;BiS2. We find a strong 
nesting of the Fermi surface at k =(tt, 7r, 0) for x=0.5 
by performing the calculation for ordered compound 
LaOo.5Fo.5BiS2- Our linear response based phonon cal- 
culation shows that the nesting results in a large phonon 
softening at this k and in the appearance of the imagi- 
nary modes associated with in-plane displacements of S 
atoms. Our \[2 x \[2 x 1 supercell total-energy calcula- 
tion finds a shallow double-well potential prompting that 
these displacements are dynamic. Contrary to the ex- 
pectations that electronic correlations may play a role in 
these systems, our calculated electron-phonon coupling 
constant A ~ 0.85 suggests that this material is a conven- 
tional superconductor. However, our anharmonic model 
calculation shows that the vicinity of the charge-density- 
wave (CDW) instability is essential for the superconduc- 
tivity which is reminiscent to ironpnictides whose prox- 
imity to SDW is well established [j. 

Our electronic structure calculations are per- 
formed within the generalized gradient approximation 
(GGA) 17]. Whenever possible, we cross-check the 
results by two different commonly used total-energy 
codes: a) the Vienna Ab-initio Simulation Package 
(VASP}|il and b) the Quantum ESPRESSO package 
(QE)[l9(. The consistency of our results for two sets of 
calculations is satisfactory. We use a 500 eV plane-wave 
cutoff and a dense 18 x 18 x 6 k-point mesh in the 
irreducible Brillouin zone (IBZ) for self-consistent 
calculations. For structural optimization, the positions 
of ions were relaxed towards equilibrium until the 
Hellman-Feynman forces became less than 2 meV/A. 
For the phonon calculations, we adopt a scalar rela- 
tivistic version of density-functional linear-response 
method [13] as implemented in QE[3|. For consistency, 
the results presented in this paper are obtained by scalar 
relativistic version of QE, unless otherwise specified. 
The effect of spin-orbit coupling (SOC) is discussed in 
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the Supplemental Material. 

Powder x-ray diffraction (XRD) pattern shows that 
LaOi_ a; F a ;BiS2 forms a layered crystal structure with 
a space group P4/nmm[l|. La, S and Bi locate at 
2b position, while O/F take the 2a site. Similar to 
LaOi-^F^FeAsfJ], the structure consists of alternating 
La(Oi_ x F x ) and BiS2 layersfH- One BiS2 layer contains 
two BiS planes (namely Bi-Sl plane) and two pure S 
planes (i.e. S2 plane) as shown in Fig. la. To study 
the influence of F doping, we carried out calculations 
for two systems, i.e. LaOBiS2 (x=0) and LaOo.5F .5BiS2 
(x=0.5). Being embedded into LaO plane, we expect that 
the substitution by F has only small effect on the BiS2 
layer. Thus we simulate LaOo.5Fo.sBiS2 by replacing half 
of the Oxygen 2a-sites by F orderly, despite the substitu- 
tion may be random in reality. The optimized lattice pa- 
rameters and Wyckoff positions for each atom are shown 
in Table I, together with available experimental data [l|. 
While the overall agreement between numerical and ex- 
perimental structures is good, Table I reveals an interest- 
ing aspect: the differences between the experimental and 
theoretical values of the ^-coordinate of SI are unusually 
large, similar as in the FeAs-superconductors 2lj. There 
is also a large difference between the numerical and ex- 
perimental inter-layer distance. Note that we have also 
performed the internal atomic coordinates optimization 
based on the experimental lattice parameters (shown in 
the Table IS of Supplementary Information). As can be 
seen from Table IS of Supplementary Information, the 
lattice parameters, namely the inter-layer distance has 
only small effect on the z-coordinate of SI. In Table I 
we also list the numerical data for LaOBiS2- Compar- 
ing with the results for LaOo.5Bio.5S2, one can conclude 
that F substitution has only little effect on the lattice 
parameters, and a small effect on the BiS2 layer. 

The discrepancy between the experimental and theo- 
retical positions of SI atoms results in a considerable dif- 
ference for valence band as shown in Fig. 2a. We also per- 
form calculation based on the numerical internal atomic 
coordinates and experimental lattice parameters. Com- 




FIG. 4: Calculated phonon dispersions for LaOo.5Fo.sBiS2 
using density functional linear response approach. 



FIG. 3: Calculated Fermi surface of La(Oo.sFo.s)BiS2: (a) 
cross section for k z —Q, (b) 3D view. 



paring with the results based on theoretical structure 
shows that changing inter-layer distance only slightly af- 
fects the band structure, and the considerable difference 
shown in Fig. 2a is mainly due to z-coordinate of SI. For 
LaOBiS2 both S 3p and O 2p states appear mainly be- 
tween -4.0 and 0.0 eV. Although located primarily above 
the Fermi level, Bi 6p have also a considerable contribu- 
tion to the states between -4.0 and 0.0 eV, indicating a 
strong hybridization between Bi 6p and S 3p states. In 
agreement with previous calculation[13| our results show 
that LaOBiS2 is an insulator with a band gap of 0.82 
eV. A little dispersion along T to Z line clearly shows 
a two dimensional character of the band structure which 
indicates that the interlayer hybridization is small. Com- 
paring Fig. 2a and Fig. 2b, it is clear that the main influ- 
ence of F substitution is a carrier doping characterized 
by the associated upshift of the Fermi level towards the 
Bi 6p band, and the system becomes metallic as shown 
in Fig. 2a. It is interesting to notice that doping by F has 
a negligible effect on the lowest conduction band. Thus 
we conclude that rigid band approximation is valid for 
LaOi-a-FxBiSa. 

There are four bands that cross the Fermi level and 
result in a large two dimensional-like Fermi surface that 
is shown in Fig. 3. Consistent with the tight -binding 



result 13], our density functional calculation shows a 
strong Fermi surface nesting at wavevectors near k = 
(7r,7r,0). 

The Bi 6p orbitals are spatially extended and strongly 
hybridized with S 2>p states near the Fermi energy. Thus 
we do not expect electronic correlations to be essential 
for this compound. To check whether the conventional 
electron-phonon mechanism can be responsible for su- 
perconductivity here, we first perform a linear response 
phonon calculation [2CJ as implemented in QE[19|. An 
18 x 18x6 grid was used for the integration over IBZ. Our 
calculated phonon spectrum along major high symmetry 
lines of the Brillouin zone is shown in Fig. 4 where the 
phonon dispersions are seen to extend up to 400 cm^ 1 . 
The phonon modes have only a little dispersion along T-Z 
direction, which again indicates the smallness of the in- 
terlayer coupling. There are basically three panels in the 
phonon spectrum that are easily distinguished along the 
T-Tj direction. The top four branches above 300 cm -1 are 
mainly contributed by O and F, while the branches below 
80 cm -1 come from the BiS2 layer. The phonon vibra- 
tions within the xy-plane show a significant dispersion as 
shown in Fig. 4. Analyzing the evolution of the phonon 
eigenvectors in the Brillouin zone reveals that there is 
clear separation between the xy and z polarized vibra- 
tions, and most of the modes show a definite in-plane or 
out-of-plane character. 

A striking feature of this phonon spectrum is the pres- 
ence of phonon softening around the M point that we 
associate with the strong Fermi surface nesting. We find 
that there are four totally unstable modes, mainly con- 
tributed by the SI in-plane vibrations, where they either 
displace along x or y direction, and either in-phase or 
out-of-phase between the two BiS planes. The polariza- 
tion vectors are shown in Fig. lb. 

To understand whether the CDW instability is present 
in this material, we perform a frozen phonon analysis 
with the unit cell doubled (V2 x \/2 x 1) according to 
the (7r,7r,0) nesting wave vector. We perform four calcu- 
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FIG. 5: Calculated double well potential for the unstable 
phonon mode using the frozen phonon method and its first 
three eigenstates. 
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FIG. 6: Calculated band structure. Black line is the undis- 
torted bands, red line is a distorted bands corresponding to 
the SI in-plane displacement. See text for details. 



lations by moving the atoms according to the eigenvec- 
tors of the four unstable phonon modes at the M point. 
Our frozen phonon calculations show a shallow double 
well potential where the SI atoms shift about 0.18 A 
away from the original high symmetry position as we 
show in Fig. 5. To illustrate the crucial change in the 
electronic structure due to M point frozen phonon SI in- 
plane motion, we show in Fig. 6 the band structures 
of distorted and undistorted structures in the vicinity of 
the Fermi level. The black lines depict the undistorted 
energy bands while the red lines correspond to the frozen- 
phonon distorted structure with SI displacement by 0.18 
A, For comparison, both lines are drawn in an y/2 x y2 x 1 
supercell. A considerable difference in band structures in- 
duced by the in-plane SI displacement indicates a large 
electron-phonon coupling from this CDW instability. 

The depth of the double well is only about 100 cm -1 
indicates that the displacements are dynamic. By solv- 
ing a corresponding anharmonic oscillator problem we 
deduce a ground state atomic wave function which is in- 



deed centered at the high symmetry position as shown in 
Fig. 5 together with the first three energy levels. It is 
therefore clear that our instable modes are not related to 
a statically distorted structure of LaOo.sFo 5 BiS2- This is 
consistent with experimental observations that the resis- 
tivity changes smoothly from 300 K to about 10 K where 
the superconductivity occurs [li [l^. 

We finally turn our discussion to the wave-vector 
(q) and mode (y) dependent electron-phonon coupling 
A„(q). At first, this can be done for all stable phonons. 
Calculation shows that the O/F modes have negligible 
contribution to electron-phonon coupling. With strong 
hybridization, the coupling, however, is relatively strong 
for the BiS based modes. For example, we can find 
A's of the order of 1 for the S2 based optical phonons 
around 310 cm -1 near the V point. The analysis of the 
polarization vectors shows that these vibrations involve 
S2 movements towards Bi atoms. Unfortunately, finding 
the integral value of A is a challenging problem due to 
the appearance of the imaginary frequencies, although 
neglecting completely the unstable modes results in al- 
ready large average coupling constant (0.75) calculated 
using 4x4x2 q-mesh. This is mainly due to the discussed 
S1/S2 vibrations. 

To find the contribution for the four anharmonic 
modes, we follow the strategy of Ref. (25[, where tran- 
sitions from the ground to all phonon excited states of 
the anharmonic well need to be taken into account. The 
detailed theory using total e nerg y frozen phonon method 
has been elaborated in Ref. [26[ . Our numerical value of 



A for the four anharmonic modes at the M point is 0.4. 
It has to be weighted somewhat by the area of the Bril- 
louin Zone where the actual instability occurs therefore 
adding it the result for harmonic A should give us a total 
coupling constant of 0.85. Inserting this value into the 
McMillan formula for T c , with the Coulomb parameter 
H* ~ 0.1 and u> D = 260 K yields values of T c ~ 11.3 K 
in reasonable agreement with experiment. 

Finally, we have also performed calculations to check 
the effect of SOC, and show the results in the Supple- 
mental Material. As can be seen from the Supplemental 
Material, the optimization including SOC gives very simi- 
lar lattice parameters and internal atomic coordinates as 
those without SOC. We find that SOC does affect the 
band structure a little, but it does not change the Fermi 
surface. As shown in the Fig.S2 of Supplemental Mate- 
rial, the Fermi surface nesting does not disappear in the 
presence of SOC. Moreover, the CDW instability also re- 
mains to be active, as shown in the Fig. S3 of Supplemen- 
tal Material. The band calculation with SOC also shows 
that SI in-plane motion has a large effect on the band 
structure, same as the case without SOC (Fig. S4 of Sup- 
plemental Material), indicating strong electron-phonon 
coupling. The large electron-phonon coupling can be also 
found from the comparison of band structure from theo- 
retical structure and experimental structure (Fig. 5S of 
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Supplemental Material) . 

In conclusion, electronic structure, lattice dynamics 
and electron phonon interaction of the newly found su- 
perconductor LaOo.5Fo.5BiS2 have been investigated us- 
ing density functional theory and linear response ap- 
proach. A strong Fermi surface nesting at (7r,7r,0) re- 
sults in large phonon softening and strongly anharmonic 
double well behavior of the total energy as a function of 
the in-plane S displacements. A large electron-phonon 
coupling constant A=0.85 is predicted, which empha- 
sizes that LaOo.5Fo.5BiS2 is a strongly coupled electron- 
phonon superconductor. 
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SUPPLEMENTARY MATERIAL 

We have optimized all independent internal atomic co- 
ordinates based on the experimental lattice parameter, 
and show the obtained results in Table IS. Note here for 
consistency, the results shown in Supplementary Mate- 
rial are all obtained using Vienna Ab-initio Simulation 
Package (VASP). As can be seen from the Table IS, the 
inter-layer distance has only small effect on the discrep- 
ancy between theoretical and experimental z-coordinates 
of SI atom. 

To study the influence of spin-orbit coupling (SOC), 
we have also performed the calculation including SOC. 
The optimized lattice parameters and Wyckoff positions 
of La(Oo.5F .5)BiS2 are shown in Table 2S. For compar- 
ison, we also list the corresponding values from the cal- 
culation without SOC in Table 2S. It is found that the 
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TABLE III: Lattice parameters and Wyckoff positions of 
La(Oo.sFo.5)BiS2 from the calculation with and without SOC. 



Site without SOC with SOC 



a(A) 




A n*7QQ 

4.U/82 


A H7Q 1 
4.U (Ol 


c(A) 




13.4835 


13.4566 
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La (26) 


0.1055 


0.1051 


z 


Bi (26) 


0.6161 


0.6163 


z 


SI (26) 


0.3891 


0.3831 


z 


S2 (26) 


0.8154 


0.8149 


z 


(2a) 


0.0000 


0.0000 


z 


F (2a) 


0.0000 


0.0000 



SOC has only very small effect on the lattice structure 
of La(Oo.5F .5)BiS2. Comparing with the results with- 
out SOC from QE (Table I in the manuscript) shows the 
consistency of VASP and QE is satisfactory. 

We show the band structure with and without the SOC 
in Fig. IS. As can be seen from Fig. IS, the SOC has a 
considerable effect on the bands around X point. To be 
specific, there are four bands crossing the Fermi level in 
the band structure without SOC, whereas including SOC 
pushes two bands above the Fermi level. Nevertheless, it 
is very interesting to notice that the Fermi surface from 
SOC calculation is still quite similar to that without SOC 
(see Fig. 2S). As shown in Fig. 2S, including SOC does 
not change the two dimensional-like feature of the Fermi 
surface. Moreover, the strong Fermi surface nesting can 
still be found at wavevectors near k —(tt, it, 0), as shown 
in Fig. 2S(b), and the main effect of SOC is to eliminate 
a small pocket around X point. 

To check the effect of SOC on the CDW instability, we 
also perform four SOC calculations by moving the atoms 
according to the eigenvectors of the four unstable phonon 
modes at the M point. Our frozen phonon calculations 
show that the double well potential remains, though the 
well-depth drops considerably, as shown in Fig. 3S. 

We show in Fig. 4S the band structures of distorted 
and undistorted structures due to M point frozen phonon 
SI in-plane motion in the vicinity of the Fermi level. 
Same with the calculation without SOC, considerable dif- 
ference in band structures induced by the in-plane SI dis- 
placement which again indicates a large electron-phonon 
coupling from this CDW instability. The large electron- 
phonon coupling can be also found from the comparison 
of band structure from theoretical structure and experi- 
mental structure as shown in Fig. 5S. 

As a final comment, we recall that the calculation with- 
out SOC is found for heavy element Pb to be satisfactory 
therefore we do not foresee significant changes in our 



estimation of the electron-phonon interaction when the 
SOC effect is taken into account. 
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FIG. SI: Band structure plot 
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FIG. S2: Fermi surface of La(Oo.5Fo.s)BiS2 cross section for 
k z =0: (a) without SOC, (b) with SOC and 3D view: (c) 
without SOC, (d) with SOC. 
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FIG. S3: Double well potential for the unstable phonon mode 
using the frozen phonon method from the calculation with 
SOC. 




FIG. S4: Band structure from calculation with SOC. Black 
line is the undistorted bands, red line is a distorted bands 
corresponding to the SI in-plane displacement. 



FIG. S5: Band structure from the calculation with SOC. 



